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Abstract. Recently the study of Fay- type identities revealed some new features of the DKP 
^ | hierarchy (also known as "the coupled KP hierarchy" and "the Pfaff lattice"). Those results 

/""\ ■ are now extended to a Toda version of the DKP hierarchy (tentatively called "the Pfaff- Toda 

hierarchy"). Firstly, an auxiliary linear problem of this hierarchy is constructed. Unlike the 
ON . case of the DKP hierarchy, building blocks of the auxiliary linear problem arc difference 

operators. A set of evolution equations for dressing operators of the wave functions are also 
obtained. Secondly, a system of Fay-like identities (difference Fay identities) are derived. 
They give a generating functional expression of auxiliary linear equations. Thirdly, these 
difference Fay identities have well defined dispersionless limit (dispersionless Hirota equa- 
tions). As in the case of the DKP hierarchy, an elliptic curve is hidden in these dispersionless 
Hirota equations. This curve is a kind of spectral curve, whose defining equation is iden- 
tified with the characteristic equation of a subset of all auxiliary linear equations. The other 
auxiliary linear equations are related to quasi-classical deformations of this elliptic spectral 
' curve. 

ON [ Key words: integrable hierarchy; auxiliary linear problem; Fay-like identity; dispersionless 

^^O i limit; spectral curve; quasi-classical deformation 
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q ; 1 Introduction 

This paper is a sequel of the study on Fay-type identities of integrable hierarchies, in particular 
the DKP hierarchy pQ. The DKP hierarchy is a variant of the KP hierarchy and obtained as 
a subsystem of Jimbo and Miwa's hierarchy of the D'^ type [21 [3]. The same hierarchy was 
rediscovered later on as "the coupled KP hierarchy" [4] and "the Pfaff lattice" [SJ EJ [7] , and has 
been studied from a variety of points of view El EDI IHl H21 H31 031 H51 HBJ - The term "Pfaff" 
stems from the fact that Pfaffians play a role in many aspects of this system. The previous 
study p] revealed some new features of this relatively less known integrable hierarchy. In this 
paper, we extend those results to a Toda version of the DKP hierarchy. 

The integrable hierarchy in question is a slight modification of the system proposed by Willox 
\17\ [T8] as an extension of the Jimbo-Miwa D'^ hierarchy. We call this system, tentatively, "the 
Pfaff-Toda hierarchy" (as an abbreviation of the "Pfaff ian" or "Pfaffianized" Toda hierarchy) . 
Following the construction of Jimbo and Miwa, Willox started from a fermionic definition of the 
tau function, and derived this hierarchy in a bilinear form. The lowest level of this hierarchy 
contains a 2 + 2D (2 continuous and 2 discrete) extension 

-D x D y r(s, r, x, y) ■ r(s, r, x, y) + t(s - 1, r, x, y)r(s + 1, r, x, y) 
- t(s, r - 1, x, y)r(s, r + 1, x, y) = 
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of the usual 2 + ID Toda equation and an additional 2 + 2D equation 

D x t(s, r, x, y) ■ r(s + 1, r - 1, x, y) + D y r(s, r - 1, x, y) ■ r(s + 1, r, x, y) = 

(see the papers of Santini et al. |19j . Hu et al. [20] and Gilson and Nimmo [21] for some other 
sources of these equations). Willox further presented an auxiliary linear problem for these lowest 
equations, but extending it to the full hierarchy was an open problem. We first address this 
issue, then turn to issues of Fay-like identities and dispersionless limit. 

As we show in this paper, the Pfaff-Toda hierarchy is indeed a mixture of the DKP and Toda 
hierarchies. Firstly, we can formulate an auxiliary linear problem as a two-component system 
like that of the DKP hierarchy pQ, but building blocks therein are difference (rather than 
differential) operators as used for the Toda hierarchy. Secondly, the differential Fay identities 
of the DKP hierarchy are replaced by "difference Fay identities" analogous to those of the Toda 
hierarchy \12\ I23|. Lastly, those difference Fay identities have well defined dispersionless limit 
to the so called "dispersionless Hirota equations". These equations resemble the dispersionless 
Hirota equations of the Toda hierarchy [22 1 [2" H[25j l26|. but exhibits a more complicated structure 
parallel to the dispersionless Hirota equations of the DKP hierarchy pQ. 

Among these rich contents, a particularly remarkable outcome is the fact that an elliptic curve 
is hidden in the dispersionless Hirota equations. A similar elliptic curve was also encountered 
in the dispersionless Hirota equations of the DKP hierarchy pQ, but its true meaning remained 
to be clarified. This puzzle was partly resolved by Kodama and Pierce |27j . They interpreted 
the curve as an analogue of the "spectral curve" of the dispersionless ID Toda lattice. We 
can now give a more definite answer to this issue. Namely, these curves are defined by the 
characteristic equations of a subset of the full auxiliary linear equations, hence may be literally 
interpreted as spectral curves. The other auxiliary linear equations are related to "quasi-classical 
deformations" |28[ [29] of these curves. 

This paper is organized as follows. In Section [21 we formulate the Pfaff-Toda hierarchy as 
a bilinear equation for the tau function. This bilinear equation is actually a generating functional 
expression of an infinite number of Hirota equations. In Section [3l we present a full system of 
auxiliary linear equations that contains Willox's auxiliary linear equations. A system of evo- 
lution equations for "dressing operators" of the wave functions are also obtained. The dressing 
operators are difference operators in a direction (s-direction) of the 2D lattice; another direc- 
tion (r-direction) plays the role of a discrete time variable. Section |4] deals with the difference 
Fay identities. These Fay-like identities are derived from the bilinear equation of Section [2] by 
specializing the values of free variables. We show that they are auxiliary linear equations in 
disguise, namely, they give a generating functional expression of the auxiliary linear equations 
of Section [3] Section [5] is devoted to the issues of dispersionless limit. The dispersionless Hirota 
equations are derived from the differential Fay identities as a kind of "quasi-classical limit" . Af- 
ter rewriting these dispsersionless Hirota equations, we find an elliptic curve hidden therein, and 
identify a set of auxiliary linear equations for which the curve can be interpreted as a spectral 
curve. 

2 Bilinear equations 

The Pfaff-Toda hierarchy has two discrete variables s, r £ Z and two sets of continuous variables 
t = (ti,tz, ■ ■ •), i = (ii,?2, ■ ■ ■)• I n this section, we present this hierarchy in a bilinear form, 
which comprises various bilinear equations for the tau function r = r(s,r,t,i). In the following 
consideration, we shall frequently use shortened notations such as r(s,r) for r(s,r,t,i) to save 
spaces. 
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2.1 Bilinear equation of contour integral type 

The most fundamental bilinear equation is the equation 

Z7TZ 

+ / ^ zS +r- s '-r'-4 e c(t-t', 2 ) r(s / + 1; r ' + ! t / + u-ii fV(s - l,r - 1, t - [z~ l ]~t) 
J 2m 

+ / ^ 7 z s - r - s ' +r 'e^*-*'^ 1 )r( S / , r' + 1, t\ t' + [z])t(s, r-l,t,t- [z]) (2.1) 
J 2m 

that is understood to hold for arbitrary values of (s,r,t,i) and (s',r',t',i). This equation is 
a modification of the bilinear equation derived by Willox [171 [T8] in a fermionic construction of 
the tau function (see Section [231 below) . Note that we have used the standard notations 



z 2 z 3 



k=l 



and both hand sides of the bilinear equation are contour integrals along simple closed cycles Coo 
(for integrals on the left hand side) and Cq (for integrals on the right hand side) that encircle 
the points z = oo and z = 0. Actually, since these integrals simply extract the coefficient of 
z^ 1 from Laurent expansion at those points, we can redefine these integrals as a genuine linear 
map from Laurent series to constants: 

ciz \ „ 
— ) a n z =a_i. 
2m 

n=—oo 

As we show below, this bilinear equation is a generating functional expression of an infinite 
number of Hirota equations. 

In some cases, it is more convenient to shift s and r as s — > s + 1 and r — > r + 1. The outcome 
is the equation 

/ ^^ ,+^ ''" S ~ r ~ 2e?(^ '~ t ' ^)r(s, ' r '' *' " [Z ~ 1] ' *' )r(s + r + !> * + t 2 " 1 ] > *) 

+ I ^L z s+rs'-r>-2 e S(t-t',z) T ( sl + lr > + ^ f + fe-l], f') T ( Sj V, t - [z'Kt) 
J 2m 

J 2717 

+ f ^^~ r " s ' +^ ''e C(^ -*'' ^ " 1) r( S , , r' + 1, t' , t' + [*])r(a + 1, r, t, i - [z]). (2.2) 

By changing variables as z — > z -1 on the right hand side, this equation can be converted to 
a more symmetric form as 

f ^-y'^'-^e^'-^ris', r', t' - [z" 1 }, t')r(s + 1, r + 1, t + [s" 1 ], i) 
y 2m 

= I ^ z -s'+r'+s-r-2 e at'-t,z)r sl + ^ ^ _ [ z -l]) T f 8 r + 1, f, t + [z" 1 ]) 

y 2^ 

+ / ^- S+r+S '- r '- 2 ^ l ~ l 'Ms\ r + 1, t>, i + [z^])r(s + 1, r, t, t - [z~% 
though we shall not pursue this line further. 
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2.2 Hirota equations 

Following the standard procedure, we now introduce arbitrary constants 

a = (ai,a 2 , ■ ■ ■), a = (ai, a 2 , . . .) 
and shift the continuous variables t, t', t, i in the bilinear equation (|2.ip as 

t' —> t — a, i —yi — a, t — > t + a, t — > i + a. 
The bilinear equation thereby takes such a form as 

f d« jj8 '+r'- 8 -r e -2C(a,*) T / a / / f _ Q _ r^-ll J _ 6 w t + fl + r^-ll j + & j 
J 2m 

+ f ^-.z s+r - s '- r, - 4 e 2 ^r(s' + 1, / + 1, t - a + [z" 1 ], t - a) 
x r(s — 1, r — 1, t + a — [z" 1 ], i + a) 

= / ^l J5 - , -r'-«+r e -2e(a,*- 1 ) r / a / + ! / t _ j _ ^ _ [ Z ]) T U -i irjt + a~t + a + \z}) 
J 2m 

+ / il z «"'-^+^ e 25(a,^ 1 ) r ( s ' ) r ' + 1}t _ 0) t_ s + [ Z ]) T ( S) r -l,t + a,t + a- 

7 27TZ 

With the aid of Hirota's notations 

D tn f ■ 9 = d tn f ■ g - f ■ d tn g, D- t J ■ g = d in f ■ g - f ■ d in g, 
the product of two shifted tau functions in each term of this equation can be expressed as 

r(s',r',t — a — — a)r(s,r,t + a + [z -1 ], t + a) 

etc., where A and A denote the arrays 

A = (D tl ,D t2 ,. . . , A„, • • •), D- t = (D h ,D h ,. . . , Dz n , . . .) 
of Hirota bilinear operators, A and Di their variants 

A = (A:, ^A 2 , ■ • • , ^A„, ■ • •) , D-t = (D h , l -D- t27 . . . , ±Dt n , . 

and (a, -Dt) and (a, A) their linear combinations 

oo oo 

(a, A) = a n D tn , (a, A") = Yl ® n D ^ ' 

n=l n=l 

Let us introduce the functions h n (t), n > 0, defined by the generating function 



n=0 

The first few terms read 



t? . , ... t? 



/i (t) = l, fci(*)=*i, M*) = -r + *2, M*) = T +t i*2 + t 3, 

2 O 
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The prefactors e ±2 ^ a,z \ etc., can be thereby expanded as 

oo oo 

e ±2 5 (a, 2 ) = h n (±2a)z n , e ±2C(a,,-i) = h n {±2a)z- n . 

n=0 n=0 

Similarly, the exponential operators e^^ Dt,z 1 \ etc., can be expanded as 

oo oo 

e ±C(A,,-i) = J- h n {b t )z-\ = hn(Dt)z n . 

n=0 n=0 
The bilinear equation thus turns into the Hirota form 

oo 

h n (-2a)h n+s , +r/ _ s _ r+1 (D t )e^ + ^r(s, r) ■ r(s', r') 

n=0 

oo 

+ h n (2a)h n+s+r _ s/ _ r ^ 3 (-D t )e^ + ^T(s - 1, r - 1) • r(V + 1, r' + 1) 

n=0 

oo 

= /in(-2a)/ ln _ s , +r , +s _ r _ 1 (A-)e <a ' Df>+<a ' D * > r( S - 1, r) • r(s' + 1, /) 

n=0 

oo 

+ /i n (2a)/ ln _ s+r+ ^_ r ,_ 1 (-A-)e (a ' Dt>+(d ' D * > r( S , r - 1) • r(V, r' + 1). (2.3) 

n=0 

The last equation is still a generating functional expression, from which one can derive an 
infinite number of equations by Taylor expansion of both hand sides at a = and a = 0. For 
example, the linear part of the expansion give the equations 

(-2/i n+s / +r /_ s _ r+ i(A) + /v+r'- s -r+i(A) D tn )r(s, r) ■ t(s', r') 

+ (2h n+s+r _ s ,_ r ,_ 3 (-D t ) + h s+r _ sl _ r ,_ 3 (-D t )D tn )T(s - l,r - 1) • t(s' + l,r' + 1) 
= h- s > +r > +s - r -i(Dt)D tn T(s - 1, r) • t(s' + 1, r') 
+ /i_ s+r+s -_ r -_i(-A-)A„r(s, r - 1) ■ t(s', r' + 1) 

and 

h s / +r /_ s _ r+1 (D t )Dt n T(s,r) ■ r(s',r') 

+ h s+r _ s ,_ r ,_ 3 {-bt)D ln T{s - 1, r - 1) • r(s' + 1, r' + 1) 
= (-2h n - s ' +r ' +s - r -i(Dt) + /i_ s / +r / +s _ r _i(Z)f n )D i - n )r(s - l,r) ■ t(s' + l,r') 

+ (2h n _ s+r+s/ _ r ^ 1 (-D I ) + /i_ s+r+fl /_ r /_ 1 (-D t -)£>t n )r(s,r - 1) • r(s',r' + 1) (2.4) 

for n = 0, 1, In particular, the special case of (|2.4p where s' = s, r' = r and ?i = 1 gives the 

equation 

- Ai-Dfj'Ks, r ) ' r ( s > r ) + r ( s ~ 1) r)r(s + 1, r) - r(s, r - l)r(s, r + 1) = 0. 
Moreover, specializing (|2.3p to a = a = 0, s' = s + 1 and r' = r — 1 yields the equation 

D tl T(s, r) ■ t(s + l,r - 1) + D h r(s, r - 1) • r(s + 1, r) = 0. (2.5) 
These equations give the lowest part of the whole Hirota equations. 
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2.3 Fermionic formula of tau functions 

Solutions of these bilinear equations are given by ground state expectation values of operators 
on the Fock space of 2D complex free fermions. 

Let us recall basic constituents of the fermion system, ipj, ipj (j 6 Z) denote the Fourier 
modes of fermion fields 

oo oo 
j=—oo j=—oo 

They obey the anti-commutation relations 

bPj ^kh = s jk, bl>j , iM + = Wj , V'fc] + = o. 

|0) and (0| denote the vacuum states characterized by the annihilation conditions 

(0\ipj = for j > 0, Vj|0) = for j < 0, 
(0|V>* = forj'<0, Vj|0)=0 forj>0. 

The Fock space and its dual space are generated these vacuum states, and decomposed to 
eigenspaces of the charge operator 

oo 

Hq = '■' l Pj' l Pj : (normal ordering). 

j=-oo 

The ground states of the charge-s subspace are given by 

l> s -i •••</>() 1 0) fors>0, = f(0|^-^! fors>0, 

[tp* ■ ■ ■ ijjl^O) fors<0, [{O^-! ■ ■ ■ ip s fors<0. 

H n (n E Z) denote the Fourier modes 

oo 

H n = :ipj7p* +n : (normal ordering) 

j=-oo 

of the U(l) current 

oo 

J{z) = -4(z)r{z):= X n z- n -\ 

n=—oo 

They obey the commutation relations 

[H m ,H n ] = m5 m+n fi 

of a Heisenberg algebra. 

Solutions of the bilinear equations are now given by 

r(s,r,t,t) = (s + r\e H{t) ge~ H{i) \s-r), (2.6) 

where H (t) and H (t) are the linear combinations 

oo oo 
H(t) = Y^tnHn, H(t) = ^ i n H^ n 

n=l n=l 
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of i?n's, and g is an operator of the form 

g = exp I ^ a jk-^j^Pt : + ^2 b jk-*Pj4>k- + ^2 W^j^t' 

\ j,k j,k j,k 

Note that this operator, unlike H{t) and H(t), does not preserve charges, hence the foregoing 
expectation value can take nonzero values for r / 0. 

Let us mention that Willox's original definition |17U18j of the tau function is slightly different 
from (I2.6p . His definition reads 

f(s,r,t,t) = (s + r\e H ®e B ®ge- B ®\s-r). 

This is certainly different from our definition; for example, Hirota equations are thereby modi- 
fied. The difference is, however, minimal, because the two tau functions are connected by the 
simple relation 



r(s,r,t,t) = exp ^ nt n t n r(s,r,t,t), 



\n=l 



so that one can transfer from one definition to the other freely. 
The bilinear equation (|2.ip is a consequence of the identity 



Is — r 



dz 

— (^(^5 <8> ip*(z)g + ?p*(z)g ® rp(z)g) 

/dz 
— (giP(z) ® g?p*(z) + g?p*(z) ® gip(z)) 

satisfied by the operator g. This identity implies the equation 

/ — (s' + r' + l\e H 'tP(z)ge~ R '\s' -r')(s + r- l\e H ^* {z)g e - R \s - r) 
J 2m 

+ / — Is' + r' + l\e H 'tP*(z)ge- B '\s' -r')(s + r- l\e H ^{z)ge~ B \s 

J 2tTI 

= /— (s' + r' + l\e H 'gip(z)e- R '\s' -r' )(s + r- l\e H g^j* {z)e~ R \s - 

J 2lTl 

/dz i ~ i 
Is' + r' + l\e H ' gip*(z)e~ H '\s' - r')(s + r- l\e H ' gi>{z)e 
2m 

where the abbreviated notations 

H = H(t), H' = H(t'), H = H(i), 3' = H(i') 
are used. This equation implies the bilinear equation (|2.ip by the bosonization formulae 

(s\e H ^(z) = z^e^is - l\e H ^ g ~^\ 
(s\e H ^r{z) = z'^e-^is + 1^+^) 

and their duals 

4(z)e- B ®\s) = e- B G-W\ a + l)z' ( £&*- 1 \ 
r{z)e~ B ®\s) = e~ B ^\s - l^-e-^* -1 ). 
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2.4 Relation to DKP hierarchy 

The Pfaff-Toda hierarchy contains an infinite number of copies of the DKP hierarchy as sub- 
systems. 

Such a subsystem shows up by restricting the variables in (|2.ip as 



t' = t, s = l + r, s' = l + r', 
where I is a constant. (|2.ip then reduces to the equation 

j> ^z^'-^e^'-^ril + r\ r', t' - [z* 1 ], i)r(l + r,r,t + [z' 1 ), t) 

+ I ^L z 2r-2r'-A e at-t',z) T n + r > + X y + ^ f + r^-ll J) 

J 2tti 

x t(1 + r - 1, r - 1, t - [z~\t) = 0, 

which is substantially the bilinear equation characterizing tau functions of the DKP hierarchy. 
Thus 

T {l + r,r,t,i) = (l + 2r\e H ^ge- S ®\l) (2.7) 

turns out to be a tau function of the DKP hierarchy with respect to t. 

One can derive another family of subsystems by restricting the variables as 

s = l + l- r, s' = I - 1 - r, t' = t, 

where I is a constant. (|2.ip thereby reduces to the equation 

= I ^ z -2r'+2r-2 e Z(t'-t,z-l) T n _ ^ ^ t j> _ rjwj _ f £ + r^n 
/ 27TZ 

+ / ^ z ~2r+2r'+2 e ai-t\z-^ T{l _ 1 _ jj y + X ± £ + [z])t{1 + 1 _ ^ r _ ^ f J. _ ^ 
J 27Ti 

This is again equivalent to the bilinear equation for the DKP hierarchy. Thus 

T{l-r,r,t,t) = (l\e H{t) ge- R ®\l-2r) (2.8) 
is a tau function of the DKP hierarchy with respect to i. 

3 Auxiliary linear problem 

3.1 Wave functions and dressing operators 

To formulate an auxiliary linear problem, we now introduce the wave functions 



^i(s,r,t,t,z) = z s+r € 



+rf(t,z) T(s,r,t-[z 1 },t) 



r(s,r,t,t) 



r, t, t, z) = z ^-2 S t, Z ) ± ~ L r- It - [z^],i) 

T(s,r,t,t) 

r, t, t, z) = z~^-2 e -^ z) r(s + l,r + l,t + [z^], t) 

r{s,r,t,t) 

r, t, t, z) = z-~e-*M T(S ' r f + l) (3.1) 

T{s,r,t,t) 
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and their duals 

^>i(s,r,t,i,z) = z s ~ r d 

V 2 (s,r,t,i,z) = z s ~ r e' 

m S ,r,t,i,z) =z~ s+r e 

ms,r,t,t,z)=z- s+r e 



(t,z-^) T ( s + 1 ,r,t,t - [z]) 
r(s,r,t,t) 
2 -i) T(s,r-l,t,i- [z]) 

r(s,r,t,t) 
rCt,*- 1 ) r(s,r + l,t,i+ [z]) 

r(s,r,t,t) 
■frz- 1 ) r(s - l,r,t,t + [z]) 
r(s,r,t,t) 



(3.2) 



These wave functions are divided to two groups with respect to the aforementioned two copies of 
the DKP hierarchy. When the discrete variables (s, r) are restricted on the line s = l+r, the first 
four (|3.ip may be thought of as wave functions of the DKP hierarchy with tau function (|2.7p . 
Similarly, when (s,r) sit on the line s = I — r, the second four (|3.2p are to be identified with 
wave functions of the DKP hierarchy with tau function (|2.8I) . If the tau function is given by the 
fermionic formula (|2.6p . these wave functions, too, can be written in a fermionic form as 



i£i(s,r,t,i,z 
V 2 (s,r,t,i,z 

y* 2 {s,r,t,i,z 



and 



(s + r + l\e H ip(z)ge 


R \s- 


- r) 


(s + r \e H ge~ H \s 
(s + r — l\e H ip(z)ge~ 


— r) 
B \s- 


- r) 


(s + r\e H ge~ H \s 
(s + r + l\e H ip*(z)ge 


— r) 
~ n \s 


— r) 


{s + r\e H ge~ H \s 
(s + r — l\e H ip*(z)ge 


-r) 
~ R \s 


-r) 


(s + r \e H ge~ H \s 


-r) 




(s + r + 1 \e H gip(z)e~ 


B \s- 


-r) 


(s + r\e H ge~ H \s 
(s + r — 1 \e H gip(z)e~ 


-r) 
B \s- 


3 

-r) 


(s + r \e H ge~ H \s 
(s + r + l\e H gtp*(z)e 


-r) 
~ R \s 


3 

-r) 


{s + r e^ge -5 s 
(s + r — le fl #(z)e 


— r) 
~ n \s 


-r) 


(s + r \e H ge~ H \s 


— r) 





^i(s,r,t,t,z 



As a consequence of the bilinear equation (|2. 1 1) . these wave functions satisfy a system of 
bilinear equations. Those equations can be cast into a matrix form as 



dz 
2tH 



^> 2x2 (s',r',t',t ,z) t V* 2x2 (s,r,t,t,z) 



dz 
2tH 



^2x2(s',r',t',i',z) t ^ x2 (s,r,t,t,z) 



(3.3) 



where 



^2x2(s,r,t,t, z) 



^i{s,r,t~t,z) ^l(s,r,t,i,z) 
V 2 (s,r,t,t,z) ^* 2 (s,r,t,i,z) J ' 
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V$(s,r,t,t,z 
V* 2 (s,r,t,i,z 

*i(s,r,t,t,z 
^ 2 (s,r,t,t,z 

$>* 2 (s,r,t,i,z 



^2x2(s, r,t,t,z) = 
^2x2(s,r,t,i,z) = 

*2x2( s >?">*>M) = 

Let us now introduce the dressing operators 

oo oo 

Wi = 1 + w me- nds , Vi = "me (n+2)9s , 



Wi(s,r,t,t,z 
^ 2 (s,r,t,t,z 

%(s,r,t,t,z 

^i(s,r,t,t,z 
<S 2 (s,r,t,t,z 



n=l 

oo 

W2 = y~] w 2n e 

n=0 

oo 



n=0 



(n+2)d s 



v 2 = i + Yl 



v 2n e 



nd s 



n=l 



nd a 



W ~L = Wlne 
n=0 

oo 

W 2 = Y,W2ne nd - 



nd a 



Vl = / J vi n e 
n=0 

oo 

V 2 = J2^ n e- n8 ' 



n=0 



n=0 



where w\ n , etc., are the coefficients of Laurent expansion of the tau-quotient in the wave func- 
tions (|3.1[) and (|3.2p . namely, 



r(s,r,t-[z~%i) _ i + ^ win2 _ n ^ 

oo 



T(s,r,t,t) 
t(s - l,r- l,i- [z -1 ],t 



n=l 

ll TS OO 



n=0 



and 



r(s,r,t,t) 



T(s + l,r,t,t-[z]) 
r(s,r,t,t) 

r{s,r- l,t,i- [z]) 
r(s,r,t,i) 



n=0 



r(s + l,r + l,t + [z%i) 
r(s,r,t,t) 

T(s,r,t,i) ~ 1+ h ' 



y^ j w ln z n , 



n=0 

oo 



^ W2nZ T 
n=0 



r(s,r + l,t,t + [z\) 
r(s,r,t,i) 

T(s-l,r,t,t+[z]) 
r(s,r,t,t) 



J2vi n z n , 

n=0 

oo 

^2v 2n z n . 



n=0 



The wave function can be thereby expressed as 

V a (s,r,t,i,z) = W a z s+r e^' z \ 9 a (s,r,t,t,z) = W a z s ~ r e^ l ' z ^\ 

K(s,r,t,t,z)=V a z- s - r e"^ z \ 0>* a (s, r, t, t, z) = V a z~ s+r e^ 1 '*^. 

3.2 Algebraic relations among dressing operators 

A technical clue of the following consideration is a formula that connects wave functions and 
dressing operators. This formula is an analogue of the formula for the case where the dressing 
operators are pseudo-differential operators [3D1 ED E2] • Let us introduce a few notations. For 
a pair of difference operators of the form 

oo oo 

P= Pn(s)e nds , Q= J2 1n(s)e nd % 
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let ^f(s,z) and $>(s, z) denote the wave functions 

oo oo 
^ z)=Pz s = Pn(s)z n+S , <S>(S,Z)=QZ~ S = J2 qn(s)z- n ~ S . 

7i=— oo n=— oo 

Moreover, let P* denote the formal adjoint 

oo 

(( = — oo 

and (P) s i s the "matrix elements" 

(P) S ' S =p s - s t(s'). 

With these notations, the formula reads 

f^(s',z)$( S) z) = (Pe 9 °Q*)s's = (Qe- 9 °P*) SS ,. (3.4) 

One can derive this formula by straightforward calculations, which are rather simpler than the 
case of pseudo-differential operators [30l l3"Tj [32]. 

To illustrate the usage of this formula, we now derive a set of algebraic relations satisfied by 
the dressing operators from the bilinear equation (|3.3p specialized to t' = t and i = t. Since 
these relations contain dressing operators for two different values of r, let us indicate the (s, r) 
dependence explicitly as W(s,r), etc. 

Theorem 1. Specialization of the bilinear equation ()3.3p to t' = t and i' = i is equivalent to 
the algebraic relations 

W- a (*y)e^- r+1 ^^(*,rr + y a (ay)e (r '- r - 1)a *W>(a,r)* 

= W a {s,r')e^- r ' +l ^V p {s,rY + V a (s,r')e^ r ~ r '- 1 ^W p (s,ry (3.5) 

for a, f3 = 1, 2. 

Proof. The (1,1) component of the specialized bilinear equation reads 

/dz f dz 

— *i(s',r',z)^(s,r,z) + 6 — $t(s', r', z)*i(a, r, z) 

/dz - - f dz - 

By the key formula (|3.4p , each term of this equation can be expressed as 
i^f>i{s',r\z)n{s,r,z) = {W^s^'-^^V^s.r))^, 

f^*Wy,*)M*>r,z) = (W 1 (s,r)e^- r ' +1 ^V 1 (s,r')) ssl 

= (V 1 (s,r')e( r '- r -V d °W 1 (s,r)) s , s , 
j^f>i{s',r',z)n{s,r,z) = (Wi(s, r')e (r ~ r ' +1)as Vi(s, r)) 8 > 8 , 

f^BltfJiZ)Ms,r,z) = (W 1 (s,r)e^'- r+1 ^V 1 (s,r')) ssl 

= (V 1 (s,r')e^- r '- 1 ^W 1 (s,r)) sls . 

Thus we find that the (1,1) component of the specialized bilinear equation is equivalent to (|3.5|) 
for a = (3 = 1. The other components can be treated in the same way. ■ 
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In particular, letting r' = r in (|3.5p . we obtain a set of algebraic relations satisfied by W, V, 
W, V. We can rewrite these relations in the following matrix form, which turns out to be useful 
later on. Note that the formal adjoint of a matrix of operators is defined to be the transposed 
matrix of the formal adjoints of matrix elements as 

A B \* _ f A* C* 
C D ) ~ \ B* D* 

Corollary 1. The dressing operators satisfy the algebraic relation 

Wx V x \ * ( e 9a \ / Wi Vi V V - e& 



W 2 V 2 J V ~e~ 9s A W 2 V 2 J \ e~ 9a ^ 



or, equivalently, 



Wx Vx V / e ds \(WiV!\( -e 9 ^ , 

W 2 V 2 ) \ -e~ d ° J \ W 2 V 2 J \ e~ 9 ° J ' 1 ' 

Proof. Let us examine (|3.5p in the case where r' = r. The (1, 1) component reads 

W x e da V x * + V x e~ ds W{ = Wie a 'V{ + V^^Wf. 

Among the four terms in this relation, Wie®"V±* and V\e~~ ds W* are linear combinations of 

e~ 9a ,e~ 29a , . . ., and Vie~ a "W^ and Wie da V{ are linear combinations of e ds , e 29a , Therefore 

this relation splits into the two relations 

Vie~ 9a Wi = W ie 9a Vx , Wxe 9a V{ = Vxe~ 9a , 

which are actually equivalent. In the same way, we can derive the relations 

V 2 e~ 9a W£ = W 2 e 9a V 2 * , W 2 e 9a V 2 * = V 2 e~ 9a 

from the (2, 2) component of (|3.5p . Let us now consider the (1, 2) component, which we rewrite 

as 

Wxe 9a V 2 - Vxe~ 9a W 2 = W x e 9a V 2 - Vxe~ 9a W 2 . 

The left hand side is a sum of e 9a and a linear combination of 1, e~ 9a , . . ., and the right hand 
side is a sum of e 9a and a linear combination of e 29a , e 39a , .... Therefore both hand sides should 
be equal to e 9a , namely, 

Wxe 9a V 2 - Vxe~ 9a W 2 = Wxe 9a V 2 - V x e~ 9a W 2 = e 9a . 

By the same reasoning, we can derive the relations 

W^V^ - V 2 e-' 9a W{ = W 2 e 9a V{ - V 2 e~ ds W{ = -e~ 9a 

from the (2,1) component of (13.50 . These relations can be cast into a matrix form as (13.60 
and ([371 



(|3.6p and (|3.7p may be thought of as constraints preserved under time evolutions with respect 
to t and i. Actually, the discrete variable r, too, has to be interpreted as a time variable. Letting 
r' = r + 1 in (|3.5p . we can see how the dressing operators evolve in r. 
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Corollary 2. The dressing operators with r shifted by one are related to the unshifted dressing 
operators as 



Wi(fl,r + 1) Vx(s,r + 1) 

W 2 (s,r + 1) V 2 (s,r + 1) 

Wi(s,r + 1) Vx(s,r + 1) 

W 2 (s,r + 1) V 2 (s,r + 1) 



e 9a 
e~ 9a 

e~ 9a 







,ds 



A B 

C 

A B 

C 



Wi Vi 

w 2 v 2 

Wi Vi 

w 2 v 2 



(3.8) 



where 



t(s + 1, r 

( log ~j frc I 

T{s,r + l)J tl 



B 



t(s + l,r + 1) 



C 



t(s + 1, r + l)r(s — 1, r) 
t(s, r + l)r(s, r) 
r(a-l,r) e _ a ^ 



r(s, r) t(s, r + 1) 

Proof. When r' = r + 1, (|3,5p reads 

W Q (s,r + l)e 2ds V$ + F Q (s,r + = W a (s,r + 1)?; + F a (s,r + l) e - 29s 

These equations can be cast into a matrix form as 



W%. 



Wi(s,r + 1) V~i(s,r + 1) 
W 2 (s,r + 1) V^s.r + l) 

Wi(*,r + 1) Vi(s,r + 1) 
W" 2 (s,r + 1) F 2 (s,r + 1) 







e 29 s 



1 
-1 



W{ W% 



V", 



V* 
1 2 



(3.9) 



Noting that 



v 2 * 



Wi Vi 

w 2 v 2 



W{ 



V* 
1 2 



we can use (|3.6p and (|3.7p to rewrite this equation as 



Wi V_i 

w 2 v 2 



Wi{s,r + l) Vx(s,r + 1) 
W 2 (s,r + 1) V 2 (s,r + 1) 







Wi Vx 

w 2 v 2 



Wi(s,r + 1) Vi(a,r + 1) 
W 2 (s,r + 1) V 2 (s,r + 1) 









d a 



Wl Vi 

w 2 v 2 



By the definition of the dressing operators, the left hand side of this equation is a matrix of 
operators of the form 



e 9a +wu(s,r + 1) - w u (s + l,r) + 



*e~ 9a + »e~ 29a + 



5io(s + l,r) c g a + 
u 20 (s + l,r) 
•e~ 9s + •e~ 2,9s + • ■ 



(• denotes a function), and the right hand side take such a form as 



/ wio(a,r + l) c _ 0jt 



1, r 



^10 ( 

^ 20 (s,r + 1 
\ uho(s - l,r) 

Consequently, 



e" 9 * + 



»e 5a + «e 29s + 



»e 9s + «e 29s + 



Wi(a,r + 1) Vi(s,r + 1) 
W 2 (s,r + 1) V 2 (s,r + 1) 



Wi Vx 

w 2 v 2 



A B 
C 
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and 

-d s n \ / txtl T/l \ _1 



Wi(s,r + 1) Fi(s,r + 1) W e - * W Wi Vi \ _/^4 5 
W" 2 (s,r + 1) F 2 (s,r + 1) JV e 9 ° )\W 2 V 2 ) = \C J 1 



where 



A = e 9s + w u (s, r + 1) - u,n(s + 1, r) + w ^( s ^ + l) ds ^ 

w 10 {s - l,r) 

B = _ ^(s + lr) a% c = + 1)e _ 9s = rB 20 (s,r + l) ds 

v 20 {s + l,r) w 10 {s-l,r) 

Rewriting this result in terms of the tau functions, we obtain the formulae (|3,9p of A, B, C. ■ 
3.3 Evolution equations of dressing operators 

The dressing operators turn out to satisfy a set of evolution equations with respect to the 
continuous time variables t and i as well. To present the result, let us introduce the notations 

(^)>o = ^Pne nd % (P)>o = £>„e^, 

n>0 n>0 

(p)<o = 5>- e ^ ( p )<° = E^ en9fl 

n<0 n<0 

of truncated operators for difference operators of the form 

oo 



rids 



Theorem 2. The dressing operators satisfy the equations 

W Un + Wxe" 6 ' Vy U \ = ( A n B n \{W! ?! 

W 2 , tn + W 2 e n9 ° V 2>tn J \C n D n )\W 2 V 2 

W Un Vi, u - V ie - n9 ° \ = (A n B n \ (Wx V! 
W 2 + V Un - V 2 e~ n9 ° J { C n D n \W 2 V 2 



(3.10) 



where the subscript t n in W\ y t n , etc., stands for differentiating the coefficients of the difference 
operators by t n , 

w a ,t n — — ' V a,t n — TT. > W a t n — tt; — j Va.tr, 



dt n ' a ' ln dt n ' a > ln dt n ' a > ln dtj 

and A n , B n , C n and D n are defined as 

A n = (W ie (n+1)9s V 2 *e- 9s )> + (Vie- (n+1)& W 2 *e- as )<o, 
B n = -(Wxe^+^Vfe 9 *)^ - {V ie -^ +1 ^W^e ds )< , 
C n = (W 2 e( n+l ^V 2 *e- 9 °)> + (V 2 e^ n+1 ^W^) <0 , 
D n = -(M / " 2 e-( n+1 ) & y i *e & )>0 - (V 2 e-^ +1 ^W^e 9s )< . 

Proof. Differentiate the bilinear equation f)3.3f) by t' n and specialize the variables to t' = t, 
i = i, s' = s and r' = r. This leads to the equation 

f ^^* 2x2 ( s > r '*) ' ^2x2(5^,^) = j> ^-9^*2^(8, r,z) • t ^ x2 (s, r,z). 
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The i n -derivatives of the wave functions in this equation can be expressed as 
d tn V a ( S ,r,z) = (W aU +W a e nd °)z s+r e^ z \ 
d tn K(^r,z) = (V a>tn - V a e- n9 °)z- s - r e-^ z \ 
d t J a (s,r,z)=W a , tn z s - r e^ z - 1 \ 
d t J* a (s,r,z)=V a , tn e- s+r e-tt i > z - 1 \ 



We now use the key formula (|3.4p to convert the last equation to equations for the dressing 
operators. Those equations can be cast into a matrix form as 

W 1)tn + Wxe" 8 ' V 1>tn \ ( -e 9 ° \ f W? W£ 
W 2 , tn + W 2 e nd ° V 1:tn J V e~ 9 ° )\ V? V 2 * 



W Un V Un - V ie - n9 ° \ ( -e 9 ° \ f W{ W% 
W 2 ,t n V 2 , tn - V 2 e- n9 ° J \e~ 9 ° J \ V{ V 2 * 

By ()3,6p and (|3.7p . this equation can be rewritten as 

W 1>tn + Wvf*' Vi,tn \( w i Vl n 
W 2)tn + W 2 e n9 ° V lttn )\W 2 V 2 

W lu V 1>tn - v ie - n9 ° \ ( Wt V l 
W 2U V 2)tn - V 2 e~ n9 ° )\W 2 V 2 

Let A n , B n , C n , D n denote the matrix elements of both hand sides of this equation. To complete 
the proof, we have to show that these difference operators do have the form stated in the theorem. 
To this end, let us compare the two different expressions 

A n B n \ _ ( Wxu + w ^ nds Vi, tn \ ( W x V^' 1 



and 



v W 2 , tn + W 2 e n9 ° V 2)tn J \W 2 V 2 

Wx^ + Wxe™ 9 * Vx, tn \( e 9 °Vfe- 9 ° -<?'V?e a ' 
W 2 , tn +W 2 e n9 ° ^2,t„ j \ -e- a 4 2 *e- & <r Q 'W{£* 

A n B n \ = f W 1>tn V ljtn - V ie - n9 ° \ ( Wx Vi 
C n D n J \ W 2)tn V 2jtn - V 2 e- n9 ° )\W 2 V 2 

Wx,t n Vx, tn - Vxe- n9 ° \f e 9 °V 2 *e- 9 ° -e 9 *V{e 9 ° 
W 2 , tn V 2<tn -V 2 e~ nd * J V -e- & W 2 *e- & e~ a 'Wt^' 

of the matrix of these operators that can be derived from the foregoing construction and the 
algebraic relations (|3.6p and (|3.7p . As regards A n , this implies that 

A n = (Wx,t n + Wxe n9a )e 9a VZer 9s - Vx^e^'W^e' 9 * 
= Wx, tn e 9 °V 2 *e- 9 ° - {Vx, tn - Vxe^e^W^. 
From the ( )>o part of the first line, we have the identity 

(A n )> = (^ie( n+1 ^y 2 *e- & )> , 
and from the ( )<o part of the second line, similarly, 

(A n ) <0 = (Vxe-( n+ V a °W 2 *e- 9 *) <0 . 
Thus A n turns out to be given by the sum of these operators as 

A n = (Wxe {n+1)9s V 2 *e- 9s )> + (V ie ~ {n+1)9s W^e^)^. 
The other operators B n , C n , D n , too, can be identified in a fully parallel manner. ■ 
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In much the same way, the following evolution equations in t can be derived. 

Theorem 3. The dressing operators satisfy the equations 

W h t n V Un - V n e n9 ° \ = (A n B n \fW 1 V l 
W 2 , in V 2 ,i n ~ V2e n9s J \C n D n J V W 2 % 
W liin + W ie - nd ° V lfn \_(A n B n \(W l V l 



where 



Wot + W 2 e~ nd ° V 2in ) \C n D n )\w 2 V 2 



A n = (^ 1 e-( n - 1 ^y; e - & ) <0 + (V^-^W^e- 9 ^, 
B n = -(Wie-( n - 1)9 'V?e d ')< - (V ie {n - 1)ds W?e d °) >0 , 
C n = (Ty 2 e-( n - 1 ) & y 2 *e- & ) <0 + (V 2 e^-^ a 'W^e- a -)>o, 
D n = -(W2e~ (n_1)9s ^e as )<0 - (^e (n ~ 1)9s W^)>o- 

3.4 Auxiliary linear equations 

The evolution equations (|3.10p . (|3.1ip and (|3.8p for the dressing operators can be readily cast 
into auxiliary linear equations for the wave functions. 

Corollary 3. The wave functions satisfy the following linear equations: 

8 ( *i *i n\ = f a b\(*i \ 

' % t>2 ^2 / V C 0/1*2 *2 *2 *2 / ' 

\&1 vp* i|r 1 ^>\ \ _ ( A n B n \ ( *i 



dt " ( ^ 2 $ 2 $;J V C n D n J \ ^ 2 *2 ^2 *5 ' ' 

*! *; * ! *; \ = / i n b„ \ / *; 

*n I iTt. iT(* lTf_ iT(* I I n n )Tl. iT(* lTf_ iT(* I ' l* - 1 ^ 



5r 2 f 2 \ C n D n J \ ^2 *2 *2 *2 

Note that each of (|3.13p . (I3.14|) and (I3.12|) is a collective expression of four sets of linear 
equations, namely, 

e^i \ = (AB\($ 1 \ f a t „$i \ = ( A n B n \ ( *! 

Otn^l \ = ( A n B n 

J V <?n A. J\ $2 

for the four pairs $ a = \^ Q , 'J/* , ^* (a = 1,2) of wave functions. The lowest (n = 1) 
equations of (13. 13h and (|3.14p agree with Willox's result [El [18] : 

d tl $i(s,r,/x) 
d tl $ 2 (s,r,iJ,) 

" 9 " I flQC T(s , +1 ' r) > l r( a +l,r+l) fl. \ 

r(Ar-l) e -q. _ e - ds + A r(.-l,r)\ I ^ $ 2 ( S , r , M ) 

r(s,r) V 6 r(a,r) / tl / 

%$i(s,r,/i) 
d h ^2(s,r,n) 



T(s+l,r)r{s-l,r) Q 8 t(s+1,t-)t(s r+1) Q s \ , . 

e e | f $ l( s > r >M) 

T(s-l,f)T(s,f-l) g, T(s + l,r)r(s-l,r) Q B I I $ 2 ( s r n) 



(3.16) 
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Thus we have obtained auxiliary linear equations for the Pfaff-Toda hierarchy. Apart from 
the fact that difference operators play a central role, this auxiliary linear problem resembles 
that of the DKP hierarchy (see Appendix). This is a manifestation of the common Lie algebraic 
structure [21 13] that underlies these hierarchies. 

Let us specify an algebraic structure in the building blocks of the auxiliary linear problem. 
Let U, U, P n , P n and J denote the matrix operators 

JJ _(W 1 V,\ fl _(W 1 Vi \ (4 B n 

V w 2 v 2 ) ' u ~\w 2 V 2 ^ n ~\C n D n 

p n = (i" f-V j=( ° j-* = -j. 



C n D n J ' V ~ e ~ da 

With these notations, (13.61) can be rewritten as 



u* = ju^j- 1 , u* = ju~ x j- x . 

This exhibits a Lie group structure (now realized in terms of difference operators). Moreover, 
(|3.10p and p. lip imply that P n and P n can be expressed as 

-e nds / tn \ 



We can confirm by straightforward calculations that P n and P n satisfy the algebraic relations 
p* — — jp j- 1 p* — — jp j- 1 

which are obviously a Lie algebraic version of the foregoing constraints for U and U. In compo- 
nents, these relations read 

A* — —pd" n p-d<> R* — p-®s R P -9s rt* _ p d a n p d s _ _ p -d s a p d s 

and 

A* — _pd s f) -ds p>* _ „-d a o -ds r<* _ p dsn p d s f)* _ _ p ~d s A p d s 

These relations are parallel to algebraic relations in the case of the DKP hierarchy (see Ap- 
pendix). 

4 Difference Fay identities 

4.1 How to derive difference Fay identities 

We now derive six Fay-like identities with parameters A and (i from the bilinear equation (|2.2p 
by specializing the free variables therein as follows: 



la) t' 


= t+[\- 


-i} + [^),t' = 


t, s' 


= s + 1, r' = r; 


lb) t' 


= t+[X 


- 1 ] + *' = 


i, s' 


= s, r' = r + 1; 


2a) t' 


= t,i' = 


i + [A] + H s> 


= s 


— 1, r' = r; 


2b) i! 


= t,i' = 


t+[A] + M, s' 


= s, 


r' = r + 1; 


3a) t' 


= t+[X 




s' = 


s, r' = r; 


3b) t' 


= t+[\- 


i' = i + [n], 


s' = 


s — 1, r' = r + l 



To clarify the meaning of calculations, we now assume that the integrals in (|2.2p are contour in- 
tegrals along simple closed curves and Co encircling the points z = oo and z = respectively. 
A and // are understood to be in a particular position specified below. 
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la) and lb). A and fi are assumed to sit on the far side (closer to z = oo) of the contour C 
The exponential factors in the integrand thereby become rational functions as 

Mt'-t,z) = I Mi'-lz- 1 ) 1 

(i- z /\)(i-z/ny 



Thus the bilinear equation (12.2f) reduce to the equations 



/ T~-< rv? V r ( s + W + [A" 1 ] + [/z" 1 ] - [z~\ 1)t(s + l,r + 1, t + [z~\ t) 

J Coo 2m z[z - \)[z - n) 



i — zr( S + 2,r,t + [A" 1 ] + ^- 1 ],i- [z])r(s, r + 1, t, t + [z]) 
/Co 2 ™ 

+ $ — ^--r(s + l,r + l,t + [A" 1 ] + [^ _1 ],t + [z])t(s + l,r,t,t- [z]) 

JCa 2m Z 



'C 

in the case of la) and 



/ TT? V r ( s ' r + t + t A_1 ] + - [ z_1 ]' *) r ( s + l,r + l,t+ [z- 1 ], t) 

/ Coo 27rtz{z- X){z-fi) 

+ £ ^ (Z "^~ M) t ( 8 + hr + 2,t + [A" 1 ] + [^] + [z-i], t)r(a, r, t - [z' 1 }, t) 



f d^ 1 

= * -^-r( S + l,r + l,i+[A- 1 ] + [ Ai - 1 ],t-[z])r(s,r+l,t,t+[z]) 
f dz 

+ * — zt(s, r + 2, t + [A" 1 ] + [A*" 1 ],* + [z])t(s + l,r,t,t- [z]) 
JCo 2m 

in the case of lb). The contour integrals in these equation can be calculated by residue calculus. 
For example, the first integral on the left hand side is given by the sum of residues of the 
integrand at z = A, //; the other contour integrals can be treated in the same way. The outcome 
are the equations 



A — fx 



t(s + 1, r, t + \pT% t)r{s + 1, r + 1, t + [A ],i) 

r(a + 1, r, t + [A" 1 ], t)r(s + 1, r + 1, t + [/i" 1 ], t) 



fj, — X 

+ ^r(s + 2, r + 1, t + [A^ 1 ] + [/i" 1 ], t)r(a, r, t, t) 
-ii i r.,-i 



= r(s + 1, r + 1, i + [A ] + \pT% t)r(s + 1, r, t, i), (4.1) 
-r(a, r + 1, t + [/i -1 ], t)r(s + 1, r + 1, * + [A" 1 ], t) 



A — /i 

' -r(a, r + 1, t + [A" 1 ], t)r(s + 1, r + 1, t + [/x" 1 ], t) 



/i — A 

+ T~ T ( S + 1) r + 2, t + [A -1 ] + ti-yXs.r.M) 

= t(s + 1, r + 1, t + [A" 1 ] + [/T 1 ], 1)t(s, r + 1, t, i). (4.2) 

2a) and 2b). A and /i are assumed to be inside Co (nearer to z = 0). The bilinear equa- 
tion (|2,2p turn into the equations 

/* dz 1 

f ^- r ( s " r ' * " * + W + + 1, r + 1, t + [z" 1 ], t) 



Pfaff-Toda Hierarchy 



19 



[' dz 1 

+ i ttA-t(s, r + l,t + [z~ 1 ],i+ [A] + [/j])r(s,r,t - [z -1 ],t) 

• y woo 

f dz z — — 

/Co 27 ™ (^-A)(«-A»J 

+ / ii ^- A )(^-^) r(s _ j r + j t t + [A] + [ ] + H)r(s + j t t _ [z]) 
J C 2m z 

in the case of 2a) and 

f dz 1 

f 77^-^(s, r + l,t- [z- 1 ], t + [A] + [h])t(s + 1, r + 1, t + [z" 1 ], t) 

• / woo 

+ f ^~^3 T ( S + hr + 2,t + [z~\ t + [A] + [/i])r(s, r, t - [z" 1 ], t) 

" Coo 

f dz z — — 

= f ^T~7 Ms + l,r + l,t,t + [\] + [fj]-[z])T(s,r + l,t,t+[z]) 

+ I dz_ (z-X)(z- m) r + ^ t - + + + + i, r> t , t - [z]) 

in the case of 2b). By residue calculus, we obtain the equations 

r(s,r + l,M+[A] + M)r(*,r,M) 
A 



A — fi 



T(s,r,t,t+[v])T(s,r + l,t,t+[\}) 



+ — ^r(s, r,t,t + [A])r(s, r + 1, t, t + [/x]) 

— A 

+ A/ir(s - 1, r + 1, t, i + [A] + [/x])r(s + 1, r, i, t), (4.3) 
r(s, r + 1, t, t + [A] + M)r(s + 1, r + 1, t, t) 
A 



A — /i 



r(« + l,r + l,t,t + M)r(s,r + l,t,t+ [A]) 



+ — ^-r r(s + l,r + l,t,t + [A])r(s,r + l,t,t + [//]) 
/i — A 

+ A / ur(s, r + 2, t, t + [A] + [/i])r(s + 1, r, t, t). (4.4) 

3a) and 3b). A and fi are assumed to be on the far side of and inside Cq respectively. 
The bilinear equation (|2.2p turn into the equations 

F dz A 
/coo 2 ™ z 2 (z-X) 

+ / ^%^T r ( s + l,r + 1, t + [A -1 ] + [z-\t + [/i])r(s, r, t - [z" 1 ], i) 



'Co 

in the case of 3a) and 

/" dz -A 
f Coo 2riz 2 (z-\) 



f dz z — — 

f 77- t(s + l,r,t + [A _1 ],t + [/x] - [z])r(s, r + 1, t, t + [z]) 

f dz z /x — 

+ f 7T~ -rtv + M + tA-y+M + ^ + l.r.M-M 

JCn 2 ™ 21 



r(fl-l,r + l,t + [A -1 ] - [z~ 1 ],t+ [/i])r(s + l,r + l,t+ [z^ 1 ],*) 
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+ <f ^%^T r ( s ' r + 2 ' * + t A_1 ] + * + MM*, r, t - [z- 1 ], t) 

Coo 

/* dz 1 

= P rT(8,r + l,t+[\-\i+\>l]-[z])T(8,r + l,t,i+[z]) 

Jco 2711 z \ z ~ N 
f dz 

+ <p —z(z - (i)t(s - 1, r + 2, t + [A -1 ], t + \jjt] + [z])t(s + l,r,t,t- \z\) 
Jco 2m 

in the case of 3b), and boil down to the equations 

A~ 1 r(s, r,t,i+ [/i])r(s + 1, r + 1, t + [A -1 ], t) 

- \- 1 t(s + l,r + M + [A" 1 ],t+ [^])r(s,r,t,t) 
= ^r(s + 1, r, t + [A _1 ],t)r(s, r + 1, t, i + [//]) 

- /xr(s, r + 1, t + [A" 1 ], t + [/j])t(s + 1, r, t, t), (4.5) 
A" 1 r(s - 1, r + 1, t, i + \jj])t(s + 1, r + 1, t + [A" 1 ],*) 

- A~ 1 r(s, r + 2, t + [A _1 ],t + [/x])r(s, r, i, t) 

= ^ _1 r(s, r + 1, t + [A _1 ],t)r(s, r + 1, t, t + [/x]) 

- ^ _1 r(s, r + l,t+ [A -1 ], t + M)r(s, r + 1, t, t). (4.6) 

Actually, these calculations are meaningful even if the contour integrals are understood to 
be genuine algebraic operators that extract the coefficient of z^ 1 from Laurent series. Thus we 
are led to the following conclusion: 

Theorem 4. The bilinear equation (|2.ip implies the Fay-like identities (|4.ip - (|4.6p . 

In the rest of this paper, (|4.ip - (|4.6p are referred to as "difference Fay identities". The 
structure of these Fay-like equations is similar to the difference Fay identities of the Toda 
hierarchy |22 [ 123 1 [T]. though the latter are three-term relations and given on a ID lattice. 

4.2 Relation to auxiliary linear problem 

We now show that the difference Fay identities are closely related to the auxiliary linear equa- 
tions. To this end, let us rewrite the identities in the language of the wave functions as follows. 

Theorem 5. The difference Fay identities (|4.1|) - (|4.6|) are equivalent to the system of the fol- 
lowing four equations: 

-D{\). 



e~°^ (., r, „) + A-*, (. + 1, r , „) - Z^'T^ ^ ^ ^ "> 

e u y A >T{s,r) /T{s,r) 

- r i ^+l,r)/r( S + l,r) fl(A) = 

r°Wr(v)/T(s + l,r + l) * 2 ^ + i,r,/x; u, 

r)c\N _ . s , „ . e fl ( A Ms - l,r)/r(s - l,r) , , , 

e D W$ 2 (s, r, M ) + A ^ 2 s - 1, r, M ) 7 -^$ 2 (a, r, M 

e -t>( l A; T ^ Sj r)/r(s, r) 

- A -i «° (A) r( a -l,r)/r( a -l,r) g(A) _ = 

e-^)d> l(s ,r,^ - d> l(s , r , M ) + A ^i^ ±± - l,r, M ) 

e- u WT\8,r)/T{8 — l,r) 

e" S(A V(s + l,r)/r(s,r) 



(4.7) 



(4.8) 



- 5 ( A )r(s,r)/r(s,r + 1) 



e- D(A) $ 2 (s + l,r,/i) = 0, (4.9) 
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nm, , . _ , . , e D ( A V(s — 1, r)/r(s, r) , , , . 

e L> WT{s,r)/T{s + l,r) 

r)(A) / _ i r )/ r ( s r ) fi . 

" D(A) w / y ()$ i g ~ W) = (4-10) 
e^(. A ir(s, r)/r(s, r — 1) 

/or i/ie /our pairs = , ^q,, ^* (a = 1,2) o/ iuai;e functions, where D(z) and D(z) 

denote the differential operators 

D(z)=f^—d tn , D{z) = f^ Z -d- tn . 

n=l n=l 

Proof. One can derive the four difference Fay identities from (|4.7p - ([4.10|) by straightforward 
calculations. Actually, this turns out to be largely redundant, namely, each difference Fay iden- 
tity appears more than once while processing the twelve equations of (|4.7p - (|4.10p . Nevertheless 
the calculations on the whole, are reversible, proving the converse simultaneously. Since the 
whole calculations are considerably lengthy, let us demonstrate it by deriving (|4.ip and (|4.2p 
from (14. 7|) and (|4.8|) for & a = ^> a (a = 1,2); the other cases are fully parallel. Recall that ^a's 
can be expressed as 

t(s, r, t, t) 

* a (s,r,z) = r(s-l,r-l,t-[z^],i) zS+r ^ {t , zl 
t[s, r, i, t) 

(|4.7p thereby reads 

e ^l A Jr(s, r)/T{s + 1, r + 1) 
Noting the identity 

e -oW e CM = (1 _ /4/A)e ((t,ri 

we thus obtain the equation 

1 _ y\ T(s,r,t- [A" 1 ] - [z^ 1 ],*) /xr(s + l,r,t - [/x" 1 ],*) 



A/ r(s,r,t - [A- 1 ],*) A r(s + l,r,t,t) 

T( 8 + l,r,t-[A- 1 ],t)T(8,r,t-^ 1 ],t) 
r(a + l,r,t,t)r(s,r,t-[A- 1 ],t) 
_ J_ ( i _ r(s + 1, r + 1, t, t)r(g, r - 1, t - [A" 1 ] - = Q 

Ap V A/ r(s + l,r,t,t)r(s,r,t- [m -1 ],*) 

for the tau function. Rewriting this equation as 

t(s, r, t - [A" 1 ] - IfT 1 ], t)r{s + 1, r, t, i) + — ^—t(s + l,r,t- i^ 1 }, i)r{s, r, t - [A" 1 ], t) 
A r..-H 



X — fJL 



T(8 + l,r,t-[\- L ],t)T(8,r,t-\pT L ],i) 



_L r ( s + i, r + i, t) t)r(s, r - 1, t - [A" 1 ] - [jr 1 ], t) = 
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and shifting the variables as t — > t + [A x ] + [// x ] and r — > r + 1, we arrive at (|4.2p . Let us now 
consider (|4.8|) . namely, 

e ^Wr(s, r)/r(s, r) 
- A -i « D(A) r(«-l,r)/r( a -l,r) D(A) _ = 
e fl Wr( S ,r)/r( a -l,r-l) U ' 

This equation turns into the equation 

A 



A — jJL 



t(s - 1, r - 1, t + [A^ 1 ] - t)r(s - 1, r, t, t) 



+ ^r(s - 2, r - 1, t - [/i- 1 ], t)r(s, r, t + [A" 1 ], t) 

- r(s - 1, r, t + [A _1 ],i)r(s - 1, r - 1, t - \pt~\i) 

M r(s - 1, r - 1, t, t)r(s - 1, r, t + [A" 1 ] - [/i" 1 ], t) = 



A — fj, 



for the tau function, and upon shifting the variables as i — > t+[fi x ], s — > s + 2 and r — ► r + 1, 
reduces to (14.11). It will be obvious that these calculations are reversible. ■ 



Expanded in powers of A, (|4.7p - (|4.10p generate an infinite set of linear equations for the wave 
functions $ a (s,r,//). As we show below, these linear equations are equivalent to the auxiliary 
linear equations (|3.13p and (|3.14p . Thus (|4.7p - (|4.10p turn out to give a generating functional 
expression of these auxiliary linear equations. 

Since the four equations (|4.7p - (|4.10p can be treated in the same manner, let us illustrate the 
calculations in the case of (|4.7p . Among the four terms in this equation, the first and fourth 
terms can be readily expanded by the identity 



e -D{X) 



°° / 1 1 \ 

J2h n (-d t )\- n , d t = id tl ,-dt 2 ,...,-d tn ,...j 



n=0 

As regards the coefficients of the second and third terms, we can rewrite them as 
e^Wr(s, r)/T(s, r) \ T(s,r) J t 

and 

e- fl W r (a + l,r)/r(a + l,r) 



e - D W T (s,r)/T(s + l,r + 1) 
r(a + l,r + l 



, 1 - log 

r(s,r) 

where / n 's denote the coefficients of the expansion 

r(s + l,: 



exp 



( (e~ D ^ - 1) log T(S . + Y ) N ) = e XP ( £ A-^ n (-^) log ■ 

oo 

= 1 + J2fn\~ n . 



n=l 
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Thus, expanding (|4.7p in powers of A, we obtain the equation 
-d tl ^( S ,r^) + (e 9 * + (log T{S T ± ^ r) ^ \ *i( a ,r,/x) 

- T(fl+ /'V 1) ee '^,r,/x) = (4.11) 

from the A -1 terms and the equations 

h n+ i(-dt)^i(s,r,fj,) - / n+ i$i(s,r,/i) 

- !li±ilL±l) £ /^^(-d^^r,//) = (4.12) 

r ^' r ^ m=0 

for n = 1,2,... from the A""" 1 terms. 

In the same way, we can decompose (|4.8p . (|4.9p and ([4.10p into and the equations 

9 tl $ 2 ( S ,r,//) + f e- a - - (log T( '~y^ j * 2 ( a ,r,M) 



r(s-l,r-l)__ a 



e- a ^ l (s,r, t i) = 0, (4.13) 



r(s,r) 

r(g + l,r)r(s - l,r)^_ a 
r(s,r) 



-a fl $i(s,r,/j) + ^ ^f-^- ds $i(s,r,/j) 



r(s + l,r)r(s,r + 1) a , . . , . -, ,\ 

-^e*$ 2 (s, r, //) = 0, (4.14) 



r(s, r) 2 

n ^ / \ t(s — l,r)r(s + l,r) fl _ . . 
r(s,r)^ 

t(s, r) z 

and the equations 

h n+ i(d t )$ 2 {s,r,(j,) - g n+1 $ 2 (r,s,ti) 

- r ( s - 1 ' r - 1 ) V 5m / in - m (5 t )e- a ^ 1 ( S ,r, / x) = 0, (4.16) 

^ +1 (-a f )$i( S ,r, M ) + T ( s + W- 1 - r ) / ne -9. $l(sr/j) 

r(s, r) z 

" T(S + 1 ' r)T( ' 2 ' r + 1) £ /m/ln-mC-^^^r.li) = 0, (4.17) 

h n+1 (d^ 2 (s,r^) + T(5 + 1, , r)r( '~ 1,r) g w e 8 ^ 2 ( g ,r, j ») 

r(s, r) z 

- r(S ~ 1 ' rM '' r ~ 1) £ ^_ m (%)e- a '$i( a> r,/i) = (4.18) 

T(s,r) z ^-^ 
m=0 

for n = 1, 2, . . ., where g n , / n and g n are the coefficients of the expansion 



ex P ((e^)-l)lo g I<^ii!i)=l + f 9 „A 
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oo 

n=l 
oo 

i+X>a" 

n=l 

Among these equations, the lowest ones (|4.1ip . (|4. 13j) . (14.141) and (|4.15p can be cast into 
a matrix form, which agrees with the lowest members (|3.15p and (|3.16p of (|3.13p and (|3.14p . 
Though the other equations (|4.12p . (14,16p . (14.17P and (|4.18p do not take such an evolutionary 
form, they can be recursively converted to the form of f)3. 13|) just as in the case of the DKP 
hierarchy pQ. Thus the auxiliary linear equations can be recovered from the difference Fay 
identities. 

5 Dispersionless limit 

5.1 Dispersionless Hirota equations 

As in the case of the KP and Toda hierarchies [33], dispersionless limit is achieved by allowing 
the tau function to depend on a small parameter (Planck constant) H and assuming the "quasi- 
classical" behavior 

T H (S, r, t, t) = e f>- 2 F(s,r,t,i)+0(H-i) {h ^ Q) (5 x) 

of the rescaled tau function 

t h (s, r, t, t) = r(/i~ 1 s, h~\ h~H, h~ l t). 

As we show below, the difference Fay identities for the rescaled tau function Tn,(s,r,t,i) turn 
into differential equations for the F-function F(s,r,t,t). Following the terminology commonly 
used in the literature, let us call those equations "dispersioness Hirota equations". 
Let us first consider (14. ip . Upon multiplying both hand sides by 

T{s + l,r+l,t,t) 

r(s + 1, r + 1, t + [A" 1 ], i)r(s + 1, r + 1, t + \fjL~ 1 ], t)r(s + 1, r, t, i) ' 

this equation turns into such a form as 

t(s + 1, r + 1, t + [A" 1 ] + [fT% i)r(s + 1, r + 1, t, t) 
t{s + 1, r + 1, t + [A" 1 ], t)r(s + 1, r + 1, t + [/i -1 ], t) 

A t(s + 1, r, t + [A~ 1 ],t)r(g + 1, r + 1, t, t) 
~~ A - n r(s + l,r + 1, t + [X~ 1 ],1)t(s + 1, r, t, t) 
fi t(s + 1, r, t + [/j- 1 ], tyrjs + 1, r + 1, t, t) 
A - n t{s + 1, r + 1, t + \pr x ] t t)r{3 + 1, r, t, t) 

1 r(s + 2,r + M+ [A' 1 ] + [/j- 1 ], t)r(g, r + 1, t, t) r(a, r, t, 1)t{s + 1, r + 1, t, t) 
+ A/i t(s + 1, r + 1, t + [A~ 1 ],t)r(s + 1, r + 1, t + [m -1 ], *) r(a, r + 1, t, t)r(a + 1, r, i, t) ' 

We can rewrite both hand sides as 

LHS = exp ( (e D W - l) (e D ^ - l) log r(s + 1, r + 1)) 

and 

RHS = — ^— exp ((e" a '' - l) (e D ^ - l) log t(s + 1, r + 1)) 



exp 



exp 



( e -SW-l)lo { 



(e B W-l)lo^ 



r(s + 1, r) 
r(s,r) 

t(s — 1, r) 



t s,r 
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- -^exp((e- 9r -l)(e D ^ -1) log r(s + 1, r + 1)) 
A — fj, 

+ J_ exp ( (e d °+ D ^ _ i) ( e 4+%) _ i) log r ( a + 1, r + 1)) 
x exp ( (e _9r - l) (e" 9s - l) log r(s + 1, r + 1)) . 



Now rescale the variables s, r, t, t as 



s -> h~ l s, r -> fT x r, t -> fT 1 *, t -> /T 1 *. (5.2) 

The last equation thereby becomes an equation for the rescaled tau function t%, in which the 
derivatives are rescaled as 

d s -f d r -f /id r , D(z) -» &D(z), S(z) -» &D(z). (5.3) 



Under the quasi-classical ansatz (|5.ip . we can take the limit of this equation as H — > 0. The 
outcome is the equation 

\ p -d r D{\)F llp -dfD{ji)F I 
e D{\)D{^)F = Ae ZJH + _L e (ds+D(\)(d a +D(ri)F+drdsF _ /g^x 

A — /i A/i 

In much the same way, we can derive the following equations from the other differential Fay 
identities (IQ)-(l4"lfl): 



e 



\ p -d s D(X)F np -d a D{ii)F i 
D(X)D(ja)F = Ae ~ + J_ e (9 r + J D(A))(9 r +D( M ))F+a r 9 s F ) /ggx 

A — /i A/x 

^d(\)d()i)f = A e w ~/i e + A ^ e{ _ 0s+D(A))( _ 9s+D(/l))F _ er9sF ^ j 

A ^ — /i ^ 

e D(\)D^)F = ~ M g + A ^ e (a r +D(A))(0 r +D( At ))F-9 r a s F ) / g 

A -1 — /i _1 

e D(\)D(^)F _ e -(d r +d a )D(ii)F _ ^ -( k a r D(x)+a B D(jj)+a r a,)F 

+ A/ze { - &+£,{A))( " 9s+B(At))F " {ar+9s)9sF , (5.8) 

e D{\)D{^)F = 1 _ V e d a {d a +D(\)-D([£))F + M e (9 r +D(A))(a r +D( At ))F_ /gg) 

A A 

The six dispersionless Hirota equations can be divided to two distinct sets { (|5.5p . (|5.7p . (|5.9p | 
and { (|5.4p . (|5.6p . (|5.8p |. The first set of equations may be thought of as analogues of the dis- 
persionless Hirota equations of the Toda hierarchy \12\ [24"1 [25| 126] . The only difference is the 
presence of the last term on the right hand side of each equation. The second set of equations 
have no counterpart in the Toda hierarchy. 

We can use {ESJ) and §EM to eliminate e D ^ D ^ F , e D ^ D ^ F and e D ^ D ^ F from 

(]5.5p . (15. Tf) and (|5.9p . This leads to the following remarkable observation. 



Theorem 6. The dispersionless Hirota equations (|5.4p - (j5.9p imply the following equations: 

^ e -dsD(\)F _ p -d r D(X)F\ + ^-1 ^ e ds(d r +d a +D{\))F _ e dr(d r +d s +D(X))F\ 

= (d r -d s )d tl F, (5.10) 
^-l^ e d a D{\)F _ e -d r D(\)F^ + ^^ e -a s (0 r -9 s +5(A))F _ e 9 r (a r -9 s +f>(A))F\ 

= (d r +d s )d h F, (5.11) 
(<9 r + c> s )c> fl F = e-^idr - d s )d tl F. (5.12) 
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Proof. We can rewrite (15.41) and (15.51) as 



/ 1 \ \ P "drD(X)F np -d r D(\)F 

e D{\)D{^)F I 1 _ _>_ e d s (d r +d s +D(\)+D(f,))F \ _ Ae ~ V e 

V V / x — n 

( 1 \ \ p -d s D(X)F _ -d s D<ji)F 

e D(X)D(ji)F I 1 _ _J_ e d r (dr+d a +D(\)+D(n))F \ = Ae V e 

Eliminating e D ^ D ^ F from these equations yields the equation 

fl _ 2_ e d s (d r +d s +D(X)+D(,j,))F\ ^ Xe -d s D(X)F _ ^g-ftflMFj 

= ^1 - -L e 9 r(dr+ds+D(X)+D( M ))F^j (\ e -d r D(\)F _ ^ e -d r D(n)F^ 

which can be expanded as 

Xe -d a D(X)F _ ^ e -dsD{tj)F + y-l e d a {d r +ds+D{\))F _ -l e d B {d r +a s +D{n))F 

= Ae~ 9rD(A)F - fie~ 9rD ^ F + X^edridr+ds+DiX^F _ -l e d r (d r +d s +D(n))F 

We can separate A-dependent and //-dependent terms to each side of the equation as 

^ e -d s D(\)F _ e -d r D(\)F^ + ^-l^ e 9 s (a r +9 s +D(A))F _ e d r (d r +d s +D(X))F\ 

= [l(e- 9aDMF - e ~ d rDMF^ + ^-lt e d s (d r +d e +D{»))F _ e d r (dr+9 3 +D(ti))F\ 

Therefore both hand sides of the last equation are independent of A and /x. Letting X, fi — > oo, 
we can readily see that this quantity is equal to (d r — d s )dt 1 F. Thus we obtain (|5.10p . In much 
the same way, we can derive (|5.1ip from (|5.6I) and (|5.7|) . If we start from (I5.8() and (15.9|) . we 
end up with the equation 

^-\^ e d s D{^F _ e -d r D{^)F^ + r-8 s {d r -d s +D{ti)F _ e d r (d r -d s +D( l j,))F\ 

_ e -d r d a F ^ e -d a D{\)F _ p -a r D(X)F\ + x -l^ e d a (dr+d s +D(\))F _ e d r (dr+ds+D(X))F^y 

By ([5TTQ]) and (joTTT]) . this equation reduces to (pU2]) . ■ 

It is easy to see from this proof that the converse is also true. Namely, if (15. 5|) . (15. 7|) and (15. 9|) 
holds, the other three dispersionless Hirota equations (|5.4p . (|5.6p and (|5.8p can be recovered 
from (|5.10p . (|5.1ip and (|5. 12[) . Thus we are led to the following conclusion. 



Corollary 4. The dispersionless Hirota equations (15.4p - (|5.9p . can be reduced to the coupled 
system of the two sets of equations <^J^, ^EM} and {(E22D, {SHIP , COD }. 



Thus the dispersionless Hirota equations of the Pfaff-Toda hierarchy can be reduced to two 
sets of equations that have quite different appearance and nature. Among the second set of 
equations, the last equation (|5. 12[) is nothing but the dispersionless limit of (|2.5p . The other 
two equations (|5.10p and (I5.11|) appear to be more mysterious. As it turns out below, they are 
related to special auxiliary linear equations for the wave functions. 

5.2 Elliptic spectral curve 

Let us examine (|5.10p . (|5.1ip and (|5.12p in more detail. Substituting A — > z, and applying (|5.12p 
to the term on the right hand side of (|5.1ip . we can rewrite (|5.10p and (|5.1ip as 

z ( e -d s D(z)F _ e -8 T D(z)F\ + z -l ^ e ds(d r +ds+D(z))F _ e d r (d r +ds+D(z))F^ _ ^ _ Q^Q^ 
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i( p a s D{z)F_ p -a r D{z)F\ + z f e -a s {dr-a 3 +D{z))F_ e a r {a r -d a +D{z))F\ = e -a r a s F^_ g^g tiF 



z [e " w — e 



Actually, these two equations are equivalent to the set of three equations (|5.10p , (|5.1ip and (|5.12p , 
because letting z —> in the second equation yields (|5.12p , and this enables one to recover (|5.1ip 
as well. 

A clue of the subsequent consideration is to introduce the auxiliary functions 

S(z) = z) + (a + r) log z - D(z)F, 

S(z) = £(t, z~ l ) + {s-r) log z + 8 S F - D(z)F. 

We can thereby rewrite these equations into the "Hamilton-Jacobi" form 

e a r s(z) = e d a s{z) + f^ Qs _ g r }Q tiF + e d s (dr+ds)F e -a a s{z) _ & d r (8 r +d s )f & ~d r S{z) ^ (5.13) 

e d r S(z) = e d a S(z) + _ Q r }Q fiF + e d s (dr+d s )F e -d a S{z) _ g 9 r (d r +d s )F & -d r S(z) ( 5 U j 

These equations imply that the exponentiated gradient vectors (e 9sS ^ z \ e drS<yZ ^) and {e daS ^ z \ 
e 8 r S(z)^ Q £ ^ e ^ wo 5_ f unc tions on the (s,r) plane satisfy the same algebraic equation 

Q = P + (0 S - d r )d tl F+ e 9s{dr+d a )F p-l _ e dr{d r +d s )F Q-l _ (5.15) 

This equation defines a curve of genus one on the (P, Q) plane. 

We now argue that this equation can be identified with the characteristic equation of a set 
of auxiliary linear equations. In other words, this curve is indeed a "spectral curve". 

The functions S{z) and S{z) may be thought of as the phase functions of the quasi-classical 
(WKB) ansatz 

*!(*) = e ft- 1 ^)+0(l) ) q, i{z) = e h-iS( Z )+0(l) {h ^ Q) 

of the wave functions ^i(z) and ^i(z). Consequently, it will be natural to expect that (15.130 
and (|5.14p are Hamilton-Jacobi equations of linear equations of the form 

e^^z) = (e w ° + o + be~ m + ce'^^z) (5.16) 

for <&i(z) = ^fi(z), ^fi(z), and that the coefficients have the quasi-classical limit 

lim a = (8s - dr)F, lim b = e 9 ^ dr+ds)F , lim c = - e 9 ^ +9 ^ F 
h-*o h-^o ft->o 

so that (|5.15p can be interpreted as the characteristic equation of (|5.16p . 

We can indeed derive such linear equations from (|3.12p . To simplify notations, let us consider 
the case where h = 1; one can readily transfer to the ^.-dependent setup of the quasi-classical 
ansatz by rescaling the variables and derivative as (|5.2p and (|5.3p . Written in terms of compo- 
nents, (|3.12p consists of the linear equations 

e dr $i(z) = A^t(z) + B$ 2 (z), e dr $ 2 (z) = C& x {z) 
for the four pairs $ a (z) = ^ a (z),^* a (z),^ a (z),^(z) (a = 1,2), where 

A = e a s+ ( log ±±lA) 1 r( 8 + l,r + l)r( 8 -l,r) ,_ 
V r(s,r + l)y tl r(s,r + l)r(s,r) 

B _ r(s + l,r + l) c g c _ r(a-l,r) c _ a , 
r(s, r) t(s, r + 1) 
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The second equation implies that $>i(z) and $2 (z) are connected by the relation 

T(s,r) 

which one can see from the definition (|3.1[) and (|3.2|) of the wave functions as well. Using this 
relation, we can eliminate $2(2) from the first equation and obtain the equation 

e^^z) = (A + Be-^C^z). 

After some algebra, this equation boils down to fl5. 16f) . Moreover, though we omit details, the 
coefficients turn out to take such a form as 

a - (\o r ( s + 1 ' r A b _ t(s + l,r + 1)t(s - l,r) 

V T(s,r + l)/ ti ' r(s,r + l)r(s,r) 

t(s + 1, r + 1)t(s, r — 1) 
t(s + 1, r)r(s, r) 

It is easy to see that these coefficients do have the anticipated quasi-classical limit. 

The other auxiliary linear equations (13.131) and (I3.14|) . too, can be converted to scalar linear 
equations of the form 

d t Mz) = (A n + Bne^C^z), dt n $i(z) = (A n + B n e- dr C)^(z) . 

Note that the 2D difference operators 

K n (e 9 ° ,e 9 r) = A n + B n e~^ C, K n (e 9 ° ,e 9 ^) = A n + B n e~ 9 ^ C 

show up on the right hand sides of these equations. Consequently, the S-functions S(z) and S(z) 
satisfy the Hamilton-Jacobi equations 

d tn S(z) = K n {e 9 ^ z \e 9 ^), d~ tn S(z) = K n (e 9 ^ z \ e 9 ^) , (5.17) 
d t J(z) = K n (e 9 ^ z \e 9 ^), d { J(z) = K n {e 9 ^ z \ e 9 ^) . (5.18) 

Those Hamilton-Jacobi equations and the constraints (|5.13p and (|5.14p may be thought of as 
defining "quasi-classical deformations" |28[ [29] of the spectral curve (15.150 . 



5.3 Comparison with DKP hierarchy 

Let us compare the foregoing results with the case of the DKP hierarchy [TJ [27] ■ 

The situation of the DKP hierarchy is similar to the KP hierarchy rather than the Toda 
hierarchy. The role of the difference Fay identities (j4.1|) ()4.6[) are played by the differential Fay 
identities 

r(r,t + [A' 1 ] + IfiL-^Tfat) 1 r(r + l,t + [A -1 ] + [^~ 1 ])r(r - 1, t) 

r(r, t + [A"i])r(r, t + ^-i]) A V r(r , t + [A"i])r(r, t + [^}) 

= 1 _ 1 dtllog f «+M) (5.19) 
X-fx T[r,t+\p !]) 

A 2 Tjr-l^ + lX-^T^t + l^ 1 }) 1? r(r,t + [A" 1 ])r(r-l,t + [ At " 1 ]) 

A - p. r(r, t + [A" 1 ] + [^~ 1 ])r(r - 1, t) A — /i r(r, t + [A" 1 ] + [ / u" 1 ])r(r - 1, t) 

= A + ,_ 8( , log lM±q±^l). (5.20, 

r(r - l,t) 
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In the dispersionless limit, they turn into the dispersionless Hirota equations 

e D{X)D(jx)F _ A -2 /i -2 e (0 r +D(A))(9 r +D( At )F = j _ 3ti (-P(A) - )-F ^ , g ^ 

A — 

>^ e -D(X)(d r+ D M )F_ _J^ e - {dr+ D(X))D M F = A + ^ _ ^ + Z)(A) + (5.22) 

A — /i A — 

for the F-function F(r, t). The first equation (|5.2ip resembles the dispersionless Hirota equation 
of the KP hierarchy. 

It is convenient to introduce the S-function 

S(z) = £(t, z) + 2r log z - D(z)F 

at this state. We can thereby rewrite (|5.2ip and (|5.22p as 



D(X)D(fi)F I ^ _ e 02F-S r S(A)-0 r 5(A t )\ = d l S W ~ gljV) 

' A — // 

D(\)D(jj,)F 

e d r s(x) _ e a r s(^ = d ti s(\) + fl^SOi) - 9r^F 



A — fj, 

and eliminate e D ^ D ^ F from these two equations. This leads to the equality 

{d tl S(X)) 2 - (d tl d r F)(d tl S(X)) - e 9rS ^ - e 9 r F - 9 ^ 

= (d tl S(fi)) 2 - (d tl d r F) (d tl S(fj,)) - e 9 ^ - e ^"^M 

hence both hand sides are independent of A and fi. Letting A, /i — > oo, we can determine this 
quantity explicitly. Thus we obtain the equation 



(d tl S(z)f - (d tl d r F)(d tl S(z)) - e 9 ^ - e 9 ' F - 9 ^ z) 

= -2d\F + l -d t2 d r F - ^{d tl d r Ff. (5.23) 



In other words, the partially exponentiated gradient vector (dt 1 S(z),e 9rS ^) of the 5-function 
satisfies the algebraic equation 

p 2 - (d tl d r F)p -Q- e 9 ' F Q- 1 = -2d 2 tl F + l -d t2 d r F - \{d tl d r F) 2 (5.24) 

of an elliptic curve on the (p, Q) plane. This curve was studied by Kodama and Pierce [27] as 
an analogue of the spectral curve of the ID Toda hierarchy. 

(I5.23|) and (I5.24p are the Hamilton-Jacobi and characteristic equations of a linear equation 
of the form 

e^^iz) = (h 2 8l + had tl + b + ce- m )<5>i(z). (5.25) 

We can derive this equation from one of auxiliary linear equations as follows. Let us again 
consider the case where H = 1. As in the case of the Pfaff-Toda hierarchy, we can eliminate 
^2 (z) and fy^i 2 ) fr° m the auxiliary linear equation (|A.7p by the relation 

^2(z) = T -^^e- 9 ^ l {z) 
r(r) 



30 



K. Takasaki 



that holds for $> a (z) = ^ a {z),^* a (z) (a = 1,2). The matrix equation (|A.7|) thereby reduces to 
a scalar equation of the form (|5.25|) for 3>i(z) = ^i(z), ^*(z) ■ The coefficients a, b, c can be 
determined explicitly as 

r(r + l)\ r(r+l)r(r-l) 



log ^— , c - 

r{r + ljtifa ~ t{t + l)t 2 _ T(r) tltl - r(r) fa 
2r(r + 1) 2r(r) 

log Z^±i2 (log T (r)) h + 2(log r(r)) Ml . (5.26) 

Rescaling the variables as (|5,2p . one can correctly recover the coefficients of f|5.23j) and (|5.24p 
in the quasi-classical limit: 

1 1 

lim a = -d u d r F, lim c = -e 9 ? F , lim b = — d t „d r F + -(d u d r F) 2 + 2d? F. 

6 Conclusion 

We have thus obtained the following equations that characterize various aspects of the Pfaff- 
Toda hierarchy: 

• the algebraic constraints (|3.6p . (|3.7p and the evolution equations (|3.8p . (|3.10p . (|3.1ip of 

the dressing operators, 

• the auxiliary linear equations (|3.12|) - (|3.14p . 

• the difference Fay identities (|4.ip - (|4.6p . 

• the generating functional expression (|4.7p - (|4.10p of the auxiliary linear equations (|3.13p 
and (BHD , 

• the dispersionless Hirota equations, (|5.4p - (|5.9p . 

• the defining equation (|5.15p of the elliptic spectral curve, 

• the Hamilton-Jacobi equations (f5T3j) . ([5TT4"]) . ([5TT7]) and (^T8|) . 

All these equations have counterparts in the DKP hierarchy. We have thus demonstrated that 
the Pfaff-Toda hierarchy is indeed a Toda version of the DKP hierarchy (or a Pfaffian version 
of the Toda hierarchy). 

Actually, this is not the end of the story. Let us note a few open problems. 

Firstly, although the other auxiliary linear equations (|3.13p and (|3. 14[) have been encoded to 
the difference Fay identities, the status of the remaining equation (|3.12|) is still obscure. Since 
its counterpart in the dispersionless limit are (|5.13p and (|5.14p , and these equations are obtained 
from the dispersionless Hirota equations, it seems likely that (I3.12p . too, can be derived from 
the difference Fay identities. Unfortunately, we have been unable to find a direct proof. If this 
conjecture is true, it leads to an important conclusion that the difference Fay identities are, 
on the whole, equivalent to the Pfaff-Toda hierarchy itself, as it is indeed the case for the KP 
hierarchy [H [33] and the Toda hierarchy [23] . 

Secondly, very little is known about special solutions of the Pfaff-Toda hierarchy. Of course 
one can freely generate solutions by the fermionic formula. Finding an interesting class of 
solutions is, however, a nontrivial problem. A possible strategy will be to seek, again, for 
analogy with the DKP hierarchy. 
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A Auxiliary linear problem of DKP hierarchy 

The DKP hierarchy has a discrete variable r and a set of continuous variables t = (ti,t 2 , ■ ■ •)• 
The Tau function r = r(r,t) satisfies the bilinear equation 



+ / ^-/^'^e^'^rir' + 1, t' + [z~ l })r{r - 1, t - [z" 1 ]) = 0. (A.l) 
This bilinear equation can be converted to the Hirota form 

oo 

^/i n (-2a)/ ln+2r ._ 2r+1 (A)e^ Dt >r(r,i) • r(r',t) 

n=0 

oo 

+ K{2a)h n+2r _ 2r ,_ z {-D t )eS a ^r{r - 1, t) ■ r(r' + 1, t) = 

n=0 

with an infinite set of arbitrary constants a = (ai,a 2 , . . .). This is a generating functional 
expression of an infinite number of Hirota equations. The differential Fay identities (|5.19p 
and (|5.20p can be derived by differentiating the bilinear equation by and specializing the 
variables as follows: 

1) t' = t+[\- 1 } + [ t x- 1 ],r' = r; 

2) t' = t+ [A" 1 ] - In' 1 ], r' = r-l. 

The bilinear equation (|A.1|) of the tau function leads to a set of bilinear equations for the 
wave functions 

< g2 (r,M)=^- 2 e«'-' T(r " 1 : t :' g "' 1) , 

r(r,t) 

¥?(r,M) = g -ar-2 e -t(M) ^ + M + M) 

r(r,t) 

V 2 (r,t,z)-z e . 
These bilinear equations can be cast into the matrix form 

/dz 
^-^2x 2 (r', t', z) ^* 2x2 (r, t, z) = 0, (A.2) 

where 

, T t z) _ (*i(r,t,z) n(r,t,z)\ /*|(r,M) *x(r,M) 

*2x2ir,t,z)-^ Mrtz) n{rAz) )i *^ r ^ z )-{ n ( r ,t,z) * 2 (r,t,z) 

The wave functions are associated with dressing operators of the form 

oo oo 

n=l n=0 

oo oo 

W 2 = ^W^d^- 2 , V 2 = 1 + Y^V2n(-d ti r n . 

n=0 n=l 
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The coefficients are determined by Laurent expansion of the tau-quotient in the wave functions 
as 



Tp -1 + 2^^ > r(r , t) -Z. 

v ' n=l v ' n=0 



oo 



-7—^ = ^W 2n Z , — — = l + ^ V2nZ 



T { r i t ) n T~{r,t) 

v ' n=0 n=l 

The wave functions are thereby expressed as 

vp a (r, t, z) = W a z 2r e^ z \ (r, t, z) = V^- 2r e -«*'*> 

for a = 1,2. 

Various equations for the dressing operators can be derived from this bilinear equation. 
A technical clue is an analogue of f|3.4j) for pseudo-differential operators [30l [31} [32] • For a pair 
of pseudo-differential operators of the form 

oo oo 

P = Pn(x)%, Q = E toW 9 *' 

n=—oo n= — oo 

let *&(x, z) and <J>(x, z) denote the wave functions 

oo oo 

V(x,z)=Pe xz = Pn(x)z n e xz , $(x,z) = Qe- xz = £ q n {x){- z) n e~ xz . 

n=— oo n=— oo 

Moreover, let P* denote the formal adjoint 

oc 

P * = E (-^)>n(^)- 

)( = — oo 

Then one has the identity 

/7 OO / i \U OO / 

*)*(*, *) = J2(PQ*)-k-i^^- = - E^n-k-i^^-, (A.3) 
fc=0 ' fc=0 

where ( )-fe-i stands for the coefficient of d~ k ~ l of a pseudo-differential operator. 
With the aid of this formula (|A.3p , one can derive the algebraic constraint 



Wx ViX ( 1 \ / Wi Fi \ 1 / -1 



W 2 V 2 J \ -I J \ W 2 V 2 J \ 1 

the discrete evolution equation 

Wx(r + 1) Vx{r + 1)\( dl \ = /4 B \( W x V 
W 2 (r + 1) V 2 (r + 1) J \ J \ C J \ W 2 V 2 

and the continuous evolution equations 

( W Un + Wid tl V Un - V 2 (-d tl ) n \ = (A n B n \(Wx Vy 

V w 2 , tn + w 2 d tl v 2U -v 2 (-d tl ) n ) \c n D n ) \w 2 v 2 



(A.4) 



(A.5) 



(A.6) 



Pfaff-Toda Hierarchy 



33 



from the bilinear equation (|A.2[) of the wave functions. Here A, B, C are differential operators 
of the form 

A = dl+ad tl +b, B = - T ±±p-, C=-^±-, 

T{r) t(t + 1) 

where a and b are the same quantities as shown in (|5.26|) . A n , B n , C n are given by 

A n = (Wxd^V 2 * + Vi(-d tl ) n WZ)> , B n = -(Wid^V? + Vi(-d tl ) n W?)> , 
C n = (W 2 d^V 2 * + V 2 (-d tl ) n W 2 *)> , D n = -(W 2 dlV{ + V 2 (-d tl ) n W?)> , 

where ( )>o stands for the projection onto nonnegative powers of dt x . These operators satisfy 
the algebraic relations 

A n = —D n , B n = B n , C n = C n , D n = —A n , 

which may be thought of as Lie algebraic counterparts of the constraint (|A.4p for the dressing 
operators. This algebraic structure is generalized by Kac and van de Leur [3] to multicomponent 
hierarchies. Let us mention that these algebraic relations among A n , B n ^ Cn: B) n are also derived 
by Kakei [9] in a inverse scattering formalism. 

The evolution equations (|A.5|) and (|A.6P can be readily converted to the evolution equations 

ttj \ _ f A B\(V 1 *5 



and 



9t «{* 2 n)-{Cn D n ){* 2 **J' (A - 8) 

for the wave functions. 
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